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A variational formulation for deriving bounds on the capture (pull-out) 
range of phase-controlled oscillators is developed. This technique is applied 
to the more common types of systems, viz., those involving symmetric com- 
parators and simple lag-RC filters. Exact asymptotic capture expressions 
relating to the simple RC filters with small bandwidth are obtained. 

I. INTRODUCTION 

The synchronization behavior of phase-controlled oscillators has been 
for many years the subject of extensive analytic studies. Although yield- 
ing significant and experimentally consistent results, these studies have 
been based primarily on laborious graphical procedures and linearizing 
approximations. In this paper we discuss and apply a variational formu- 
lation with which to analyze the synchronization phenomena of phase- 
controlled systems more directly and somewhat more generally. 

The phase-controlled oscillator is represented by the block diagram 
of Fig. 1(a). The basic function of such a system is to establish and main- 
tain, within certain limits, an approximately zero difference in frequency 
between the input and local carriers; i.e., the local oscillator is constrained 
so as to follow variations in the frequency of the input carrier. As indi- 
cated, this control is effected by applying the two carriers to a phase 
comparator and linear filter in cascade, and utilizing the resultant output 
to appropriately control the frequency of the local oscillator. The com- 
parator output in most cases is a simple function of the phase difference 
<Pi — <p between the carriers, the type of function depending to a large 
extent on the particular physical operation by which the phase difference 
is measured. For example, if the frequency spectrum of <p< — <p is base- 
band, then multiplying the two carriers and filtering out the resultant 
high-frequency component yield an output which approximates that of 
an ideal' 'cosine" comparator, viz., f(<pi — <p ) ~ cosO, — tp ). In general, 
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pig i — ( a ) Block diagram of the phase-controlled oscillator; (b) RC and lag 
filters. 

f( Vi — Vo ) i s not a baseband function ; however, the filter depicted in Fig. 
1(a) is usually a low-pass lag or RC type, serving to eliminate any high- 
frequency deviations in the control branch. 

As applied here, the term synchronization relates to the difference in 
frequency between the input and output carriers ; specifically, we consider 
the system to be in a synchronous state provided this frequency differ- 
ence is either identically or asymptotically zero, i.e., 



(cfyi/dt) - (d<p /dt) = or 



(t -> 00 ). 



Despite the fact that these conditions are realized in special cases only, 
near-synchronous states are readily established by most practical sys- 
tems. Whether an exact, near, or nonsynchronous state prevails in a given 
situation depends largely on the particular form of input dtpi/dl and the 
related initial conditions. If completely general, the analysis of such 
phenomena necessarily involves both arbitrary inputs and arbitrary 
initial conditions. Owing to the intractable nonlinear formulation en- 
countered, however, a broad treatment of synchronization does not ap- 
pear possible at present except for trivial systems. 

To establish effective criteria for evaluating system performance yet 
obviate the mathematical difficulties associated with a more detailed 
analysis, we consider special classes of input functions that in cases of 
interest characterize actual inputs and lead to exact or approximate solu- 
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tions. The most acceptable as well as convenient input is perhaps the 
step function, a single "jump" in frequency of magnitude C\(t ^ 0), 
normally applied with the phase difference initially zero and the system 
initially in a synchronous state. Such functions represent the more ex- 
treme or sudden type of input deviation, and yield in the simple lag-RC 
filter cases (cf. Fig. 1(b)) autonomous formulations to which phase 
plane descriptions and techniques are applicable. We often find it possi- 
ble to consider inputs with less restrictive initial conditions, namely those 
satisfying the asymptotic relation 

(dtpi/dt) — (o) c ) — * Ci = const (t — > oo ). 

(Here, oj,. represents the "free-running" frequency of the local oscillator 
and C 2 , an asymptotic frequency jump.*) The maximum positive values 
of Ci and C\ for which the related functions restore the system to a 
synchronous state are denned as the positive "pull-out" (capture) and 
"pull-in" ranges, respectively. As measures of performance, these defini- 
tions have been used extensively, especially in cases involving sine com- 
parators and lag-RC niters. 

One essential purpose of this paper is to formulate a variational pro- 
cedure whereby one can obtain analytically bounds on the capture range. 
The method applies specifically to the second order, autonomous, non- 
linear equation, 



<it 2 



(••I) 



which for x ~ <p< — <p„ describes phase-controlled systems with lag-RC 
filters and step function inputs. Relative to the phase plane representa- 
tion given in Fig. 2, solution trajectories of the above equation for lag- 
RC systems arc of two principal types: a nonsynchronous form running 
on indefinitely, and a synchronous or capture form satisfying the con- 
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Fig. 2 — Phase plane representation of solution trajectories. 



See list of symbols, Appendix D. 
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dition (d<Pi/dt) - (d<p /dt) — > 0. As discussed previously, the basic prob- 
lem is determining values of d for which only the latter form results. 
To find these values by phase plane methods requires fairly accurate 
graphical plots of a large number of solution trajectories. Of main con- 
cern, however, is the relationship between Ci and the resultant solution 
type, not detailed knowledge of the trajectory. Utilizing this principle, 
we show below that subintervals of d can be identified as corresponding 
to either capture or noncapture trajectories by means of "test" trajec- 
tories, prescribed functions exhibiting the essential features of the actual 
solutions. The more closely the test contours fit those of the solutions, the 
more closely the derived subintervals bound the capture and noncapture 
ranges. For obtaining optimal bounds, we normally construct a class of 
appropriate test functions and vary over that class, the procedure simu- 
lating to some extent that of linear vector-space theory where variational 
formulations based on trial eigenf unctions furnish bounds on eigenvalues. 
Although the development is restricted to the second-order, autonomous 
equation given above, straightforward generalizations enable one to deal 
with higher-order, nonautonomous equations. 

In the treatment that follows, we consider first the mathematical for- 
mulation of phase-controlled systems involving lag-RC filters, then the 
variational formulation, and finally the application of the latter to sym- 
metric comparator functions, i. e., 

/(0 = -/(-0 (« - <Pi~<Po) 

/(e) =/(« + 2r) 

/(e) ^ (0 ^ e ^ tt). 

Regarding the symmetric comparator and RC filter case, several results 
arc derived: (i) bounds on the capture range under general bandwidth 
conditions; (it) exact asymptotic formulae for the capture range under 
small and large bandwidth conditions; (Hi) the fact that for capture 
solutions, the steady-state phase error lies in the interval 

— IT ^ 6 ^ 7T. 

With respect to sawtooth comparator functions, i. e., 
/(«) - E /o(« " 2*n) 
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we determine exact capture relations for RC niters and capture bounds 
for lag filters. 

II. GENERAL FORMULATION 

As in previous investigations, 1 - 2 the mathematical representation of 
the phase-controlled system shown in Fig. 1 is based primarily on three 
physical constraints : 

(i) With regard to the control of the local oscillator, the instantaneous 
frequency co (t) of the local carrier consists of a constant, free-running 
component co c which is independent of the control signal S(t) applied to 
the local oscillator, and a component which is directly proportional to 
S(t); i. e., 

«■(*)- *5 -«. + «S(0- (1) 

In general, w c is chosen so as to lie within the frequency range of the in- 
put signal. Accordingly, we represent the input as the sum of two fre- 
quency components, co c and a time-dependent deviation Aw(t); viz., 

«*(0 = ~ = u c + Ao,(0- (2) 

(ii) Regarding synchronism, it is required that the system exhibit a 
"natural" synchronous state for which &>,-(<) = oi (t) = cc c (t ^ 0) pro- 
vided only e(0) = *>,-(0) - <p o (0) = and Au>(l) = (t ^ 0). If, after 
being initially established, this state is to persist indefinitely, then 
S(t) = (t ^ 0). The latter condition requires the comparator output 
to vanish for zero error input and the linear filter to be initially inert. 
More explicitly, it is stipulated that 

S(t) = (««) =0, t^O) (3) 

/(0) = (4) 

— ) ^ m = const < oo . (5) 

ffc/«=o v 

In addition, the usual supposition of physical realizability implies the 
requisite 

1/(0 I ^ /«• = const < co (- co < e < «,). (6) 

(tit) Relative to boundary conditions, we assume that with the system 
initially in the synchronous state outlined in (ii), there is applied at 
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t = an arbitrary MO depending only on the behavior to be analyzed: 
e(0") = 0, Aw(0 = o (< < 0). (7) 

The constitutive conditions given by (1) through (7) yield 

^ = MO - « f H(t - r)Mr)\dr (t ^ 0) (8) 

at J o 

where the filter impulse response H(t), Aw(0, and /(e) are regarded as 
generalized functions. 3 For physically realizable filters and comparators, 
the convolution in the last term of this expression is bounded in a neigh- 
borhood of the point t = 0; hence, provided Aw(0 is bounded, de/dt be- 
haves similarly in this neighborhood and 

e((T) = e(0) = e(0 + ) = 0. (9) 

Considering specifically the lag filter for which 

H(t) = £-T^ ( 8 -^)\ = -[*(0 + (<* - oi)«-*>1 (/ ^ 0) (10) 

Lo-2 \S + <Ti/ J CT2 

and differentiating (8) with respect to i, one obtains 

CTi dt af m \ <T\ dt / 



- CU-) -%+ (l - P)kgf: it > o; 



efo 



1t~~k 9 



(11) 



(12) 



where 



x(0 = -e(0, C(.t) -i/(«x) 

fl . i + * * - feY. 

7TO-2 \<XJm/ 

One notes from (1), (8), and (10) that 

I coo - co c | = a|5(0| = «I#7I ^ «/» ( H(r) dr = af m . 

Thus, the local oscillator cannot establish synchronism, i.e., w< = coo , if 
M0 = wi - co c > a/ m (/ ^ 0). Anticipating that Aco/a/ m g 1 incases 
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which involve synchronism, we define the relative deviation 12(0 by the 
relation 

0(*)-^. (13) 

Equations (11), (12) and (13) can be combined to give 

g (l f = 0(0 + -~ - C{x) - kg + (1 - m-gf. (14) 
ax a i at ax 



Also, by (4), (8), (9), (11), and (12), 

L(*) 

r<ri \dt/o+ 
x(0~) = x(0) = x(0 + ) = 0. (16) 



, n ._*.(*) _«« (15) 

tt(T\ \at/o+ k 



To render the analysis of (14) to (16) tractable, we consider only 
specific classes of £2(0- These deviations, although restricted, are chosen 
so that the associated synchronization behavior is representative of the 
more general phenomena. Accordingly, we discuss the following two 
classes: 

(i) As a characterization of sudden input changes, 12(0 is assumed to 
be a step function of magnitude £, causing either a temporary or perma- 
nent loss of synchronism. The maximum value of £ for which the system 
returns to a synchronous state is denoted by £+ ; namely, 

$+ = sup {f|f £0; g(a;)->0(«-»oo)} (17) 

where 

It (t ^ 0) 

o(0 = 

[o (t < 0) 
Similarly, for negative values of £, 

f_ = inf U|{ <0; g{x) ->0 (t -» oo)}. (18) 

We define the "relative" capture range £ c by the relation 

& = *+-*-. (19) 

(ii) As a generalization of the deviation type above, 12(0 is assumed 
initially arbitrary but asymptotically such that 

12(0 —*y = const (t — * qo ). 
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Iii a manner similar to that of (i), we define the relative "pull-in" range 
y p by the relation 

7 P = 7+ - 7- ( 20 ) 

where 

7+ = sup {7 I 7 ^ 0; g(x) -> 0,12(0 -+ 7 (J -> w )} 
7_ - inf {t It <0; g(x)-*0fl(t)-*y(t-> »)}• 

Since this latter class contains the class of (t), y P serves as a lower bound 
on & . 

The capture and pull-in range as well as the steady-state phase error 
are significant as measures of the effectiveness of the system to correct 
for input deviations. This synchronizing capability is of prime impor- 
tance; however, in the treatment here solutions x(t) are not of direct 
interest. 

III. PHASE PLANE REPRESENTATION 

The phase plane is perhaps the most natural means of describing the 
behavior of the phase-controlled system. As shown in Fig. 2, the solu- 
tion trajectories in the (g,x) plane assume two essential forms: Those 
approaching singular or critical points x on the abscissa represent the 
capture phenomenon; those running on indefinitely, the noncapture or 
nonsynchronous state. The values of .t are determined directly from 
(11) and (12). 4 ' 5 For Q(t) = f (< ^ 0), 

C(x ) = £ g(x ) = 0. (21) 

Similarly, for 0(0 -»7 (* -» =° ), 

C(x ) =7, ff(so) =0 (i-> oo). (22) 

Since in both of these instances | C(x ) | £ 1, a necessary condition for 
the presence of critical points and, consequently, capture and pull-in is 
that 

|t|£l, ItI^L (23) 

As regards the nature of these points, it is shown in Appendix A that .r 
is either a stable or saddle point according as (dC/dx) Xo is positive or 
negative. Borderline cases, e.g., those involving the condition 

(dC/dx) X0 = 

and the coalescence of stable and saddle points, are not of physical 
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interest as infinitesimal variations of the system parameters eliminate the 
borderline properties. In the event that C(x) is discontinuous, the above 
criteria are applied to a regular sequence, i.e., a sequence \C n (x)} '° of 
continuous and different iable functions for which 



C n (x) -> C(x) (n -> 
the desired conditions obtained in the limit. 



), 



IV. TEST TRAJECTORIES 



A VARIATIONAL APPROACH 



At this point we describe a method with which to obtain upper and 
lower bounds on the relative capture range. With reference to the gener- 
alized phase plane portrait of Fig. 3, the capture and noncapture trajec- 
tories are regarded as terminating on two different sets of points M and N, 
respectively, and emanating from various regions of the initial state 
domain D, e.g., D M and D N . Terminal points in N represent points at 
infinity; those in M, critical points. Saddles are included in the latter 
group because a trajectory which approaches such a point, involves an 
infinite amount of time and therefore satisfies capture conditions (17) to 
(19). Stated in the most general manner, the basic problem to be treated 
here is the determination of contours r c , namely, boundaries separating 
regions in D from which either capture or noncapture trajectories origi- 
nate. For 12(0 = £, domain D becomes the g axis and contours r c de- 
generate to points. 

As a technique for locating r c , we first construct "test trajectories," 
i.e., functions h(x), which run from D to a position somewhere between 
sets M and N. The test paths are selected so that if the solution paths 
are known to lie entirely on either side of the former, then the latter must 
terminate on points in the terminal set lying on the same side. It is then 




Fig. 3 — Generalized phase plane portrait. 
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possible to show by means of comparison theorems that for certain re- 
gions of the initial state domain positioned on one side of a test trajec- 
tory, the corresponding solution trajectories are similarly positioned. 
Consequently, the regions in D so determined are identified as either 
capture or noncapture domains forming bounds on T c . 

To make the above procedure precise, we utilize the following compari- 
son theorem: Let g(x(t)) be a trajectory branch satisfying the autono- 
mous equation 

^ = F(x,g) (xeI,teJ) 
ax 

where / denotes a closed interval of x(x a ^ x ^ x b ) and J, a branch 
interval of t defined by the inverse function t = t(x) (x e I). If for a 
test trajectory h(x), 

f.< (>) F(x,h) (x € I) 
ax 

(24) 

h(x a ) 5S (^)g(x a ) 

h(x),g(x),t(x) eC (xcIJcJ) 

then 

h(x) < (>)g(x), x a < x ^ x h . (25) 

A brief proof of (24) and (25) is outlined in Appendix B. Although nor- 
mally restricted to continuous functions, this theorem can be applied to 
generalized functions through the related regular sequences. 

With g(x a ) positioned in the initial state domain and the test tra- 
jectory constructed as a family of functions involving parameters 
Xi , • • • , X„ , i.e., h(x,\), the inequalities of (24) can usually be reduced 
to the form 

g(x a ) > (<)G(x,\) (20) 

where the inequalities are sometimes reversed. Since the G function is 
known, there are obtained from (26) regions of domain D for which the 
resulting trajectories lie either above or below those of the test functions; 
if in the sense described above, sets M and N are isolated by h(x,\). 
then these regions constitute portions of the capture and noncapture 
domains to be determined. The optimal bounds on g{x a ) for the given 
family h(x,\) are derived by extremalizing G with respect to x and X: 

> Qv = inf -sup (7(.r,X) 

g{x a ) X X • (27) 

< g L = sup -inf G(x,\) 

X x 
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The quality of these bounds depends on that of the test function in that 
the closer g(x) is approximated by h(x), the better the contours r c are 
delineated. Unlike the functions associated with Liapounoff's second 
method, 4,5 h(x) is related directly to whatever properties of the actual 
trajectory are known and is constructed so as to conform as closely as 
necessary to these properties. 

V. GENERAL SYMMETRIC COMPARATOR — RC FILTER 

In this section we apply the method of test trajectories to phase-con- 
trolled systems employing RC filters and symmetric comparators, viz., 
those defined by the relations 

C(x) = C(x + 2») (n = 0, ±1, • • • ,) (28) 

C(x) = -C(-x) (29) 

C(x) ^ (0 ^ x £ 1) (30) 

where in accordance with (3) to (5) and (6) 

C(0) = (31) 

C(x) g 1 (0 ^ x £ 1). (32) 

For the RC filter a-, -> oo (p = 1) in (14) ; therefore, 

0^ = 0+ I*-c(») -kg. (33) 

ax <s\ at 

5.1 Relative Capture Range 

With 0(0 = £, (33) reduces to 

A £ - - \k - kg - C(x)} = F(x,g). (34) 

It is noted in this expression that trajectories corresponding to inputs £ 
and — £ are symmetric about the origin of the phase plane in that the 
respective solutions g$(x) and g~i(x) are related by 

f7 £ (.r) = -g^(-x). 

Consequently, both solutions exhibit the same capture characteristics, 
and 

5+ = -€- (35) 

or 

k - 2£ + = -2£_. (36) 
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As a result, only positive input deviations need be considered. Referring 
to (15), one obtains for the initial state line 



<7(0) = jU 0. 



(37) 



In order to establish prude noncapture intervals on this line, we employ 
the following test function: 



h(x) = EmiOk - 2n) (0 £x < co) 



(38) 



where 



Mi 



(*) = 



-- f C(f) 

p Jo 



dz 



1- 



(-1 *xgl) 

(1*1 ^ 1) 



= ( C(z) dz. 
Jn 



The functional form chosen for hi(x), as well as most of those to follow, 
serves two basic aims: 

(i) In addition to having the general appearance of a noncapture 
trajectory, hi(x) isolates the set of critical points on the x axis from the 
region above the test trajectory; i.e., solution trajectories which lie 
above fa(x) are bounded away from the x axis (cf. Section IV). 

(«) With respect to factor m(x), hi(x) relates closely to the algebraic 
structure of (34). 

We first verify property (i). By (28) to (32), 



p = sup 



Thus 

and 
whence 



[ C(z)dz. 

- r dz) dz ^ i 

n J(\ 



/xi(-r) ^ 



hi(x) ^ (0 £ x £ oo ). 
The comparison theorem applied to (34), (38) and (39) insures that 
^ fa(x) < g(x) (x > 0) 



(39) 
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for all values of £ satisfying the inequalities 



£ > C(x) + kh + h ^ 
ax 



(40) 



= (l - 2^5) C(x) + M,(.t) = GCr,/?,) (x ^ 0) 

Mo)^ ? (o) = |. (41) 

Since C(0) = 0, condition (41) is contained implicitly in (40), the 
latter expression corresponding to that of (26). 

Functions g(x) as trajectories bounded away from the x axis con- 
stitute noncapture solutions; hence, any £ interval consistent with (40) 
is a noncapture interval. For obtaining rough limits, we set X 2 = 2p/c 2 . 
Inequality (40) then becomes 

£ > kh(x) (x ^ 0). 

Therefore, as indicated by (27) and (35), all values of £ which satisfy 
the relation 

£ > £i = sup khi(x) = sup km(x) = k(2p) h (42) 

lie in the noncapture interval with the extremes +£ x and — & represent- 
ing exterior bounds on the capture range. Although directly significant, 
these limits are derived primarily to provide essential information in 
more general calculations. 

For a refined treatment we construct a more complex test function 
h(x): Let g (x) denote the separatrix solution which runs from the 
initial point (g (0), 0) through the last saddle in the interval nearest 
to the point x = 1, and let £ n denote the related input deviation. We 
then define 

h 2 (x) = £ M 2 (.r - 2n) (0|.t< oo) (43) 

n=0 

where 

(mix) (\ = £ , - l ^ x < 0) 

»M = I gdx) (0 ^ z ^ .to) 
,0 (x () ^ .r g 1) 

In addition to assuming the existence and uniqueness of g (x), we sup- 
pose further that g ( ,(x) |0(0g^ .ro), for if O (.r) = k/tndx/dt < 0, 
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then x(t) decreases as time increases, and the trajectory is depicted 
receding from the point x in opposition to the definition of g (x) ; thus, 

«(a:) ^0 (-1 <x g 1) 

and 

h(x) ^ (0gs<«). (44) 

Also, 

*.(0) = flfo(0) =|. (45) 

Therefore, we have the same sufficient requirements for noncapture 
values of £ as those of (40) and (41) : 

£ < C(x) + kh + h^ = G(x,h 2 ) (46) 

A.(0) =|^ff(0) =|- (47) 

To determine sup G(x,h 2 ), we first derive several conditions relative 

to x , £ , and go . With the selection of g as a proper solution, (34), 
(46) and (47) combine to yield 

G(x,g ) = £o (0 ^ a; ^ x„). (48) 

We next let x m represent the greatest critical point value in the interval 
x < x ^ 1 ; viz., 

x m = sup {Xi | C(xi) = £o ; x < Xi S 11- 

i 

Inasmuch as x Q is required by definition to be the only saddle in the 
interval x ^ x S 1, z ro must be a stable point; then, (dC/dx) Xn > 0, 
as discussed in Section III, with the result that C(x)_ > C(x m ) for 
x m ^ x < 1. However, since C(.r) as an odd function satisfying 

C(l) =i(C(l + ) +C(D) = c(0) =° 

must either approach zero or jump discontinuously to zero as a; -» 1. 
x m cannot exist ; hence, 

C(x) 7* ft(xo <x ^1), C(xo) - £c 

and 



rf.r 
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These conditions imply that 

C(x) ^ C(xo) (x ^ x ^ 1) 
whence 

sup (7(3,0) = supC(x) = C(x ) = ft. (49) 

In addition, it is noted that because of the infinite amount of time for 
the po trajectory to reach the saddle point, £ lies within the capture 
range and outside that defined by inequality (42); i. e., 

ft £ fc = fc(2p)*. 

Consequently, 

_«up <o G(.<-,M.) = supjYl - iQ C(x) + k^(x)~\ 

= -sup fc(x)] (50) 

- 6o. 

As a result of (48), (40) and (50), the noncapture range is given by 
5 > sup G(x,h 2 ) = sup G(x,to) = go . (51) 

That values of £ less than £ are in the capture range is shown by re- 
versing the inequalities in (4(5) and (47) and invoking topological prop- 
erties of the phase plane portrait: Values of £for which 

£ < G(x,g ) = ft (0 ^ x ^ Xo) (52) 

correspond to solutions g(x) which satisfy the relation 

g(x) < g (x) (53) 

a condition indicating that for such £ the associated solution trajectories 
cither terminate on stable nodes or cross the x axis at a point x = x < 
x and proceed in the — x direction. As shown by the phase portrait of 
Fig. 4, the crossover point C in the latter case must be located so that a 
saddle S lies just to the right and a focus F just to the left, for the 
trajectory would otherwise intersect either itself or a separatrix of S. 
Using the same argument proves that g(x) must in addition spiral in 
toward F. Consequently, (52) is a capture condition and ±£ are exact 
exterior bounds on the capture range. By (23) and (36) 

ft = 2£ ^ 2 (54) 
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Fig. 4 — Phase portrait of trajectories. 



or 



€+ = & =5 I- 



(55) 



Two significant results have so far been obtained: 

(i) The possibility of the system's reaching a synchronous state in a 
region outside the first domain does not exist; i.e., the steady state phase 
error x„ must be such that 



x se t\Xi I { - C(Xi); 



!> >* 



^ Xi ^ 1 



(56) 



(») The determination of the capture range depends solely on that 
of the separatrix solution go(x). 

Among the set of stable points in (56) the appropriate x t , namely, 
the exact value of x ss , can be found by the application of the above pro- 
cedure to each separatrix spanning the intervals g{/H &/fc and 
^ x ^ .To on the initial state line and x axis, respectively. For example, 
the first stable point x t = inf {x { ) is isolated by the separatrix which 
connects the points \/k and 



= mi\x j \l = C(xj); 






dC 
d 



< 0; ^ .x-y ^ .To 



in that the solutions for all £<$(£> 0) lie below this separatrix and 
Xs. = xi . Each Xi can be handled in this manner. 

In general, g cannot be determined exactly, however, as regards 
specific comparators, graphical methods are found to yield results as 
accurate as one desires because of the relatively small distance traversed 
by the associated trajectory. The use of such graphical approximations 
for obtaining the capture range is prevalent in the literature in advance 
of any formal justification. It is important to point out that even for 
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the RC filter case, interior and exterior bounds are sometimes more 
convenient and useful than the graphical solutions. 

5.2 Asymptotic Results for Small and Large k 

It is possible in the case of vanishing A- to derive explicit relations for 
both the capture range and go(x). We first consider an exterior bound 
formulation based on the test function 



fh(x) = 2 to(x ~ 2n) (x ^ 0) 



(57) 



where 



Pa 



(*) = 



[[*-2j["c(.)«*[ 



(-1 ^ x ^ x) 



(1*1 ^ 0). 



Note that 



Mi(D = 

m(x) > (-1 < x < 1) 
h-Ax) ^ (as £ 0) 
h(0) = (2 P V 

in accordance with noncapture criteria (cf. (40) et seq.). Hence, as in 
(40), (41) and (42), exterior bounds ±fe are given by the expression 

f > & = sup G(xji 3 ) = sup G(x,fx 3 ) 

= sup 2p/,'- - 2fc 2 [ C(z)dz] (58) 

l*l*i L Jo J 

= H2p) { 
which is consistent with the initial condition of (41); viz., 

fe 8 (0) = (2p) 5 ^ ff(O) =|. 

To treat the interior bound calculation, we use 

h A (x) = [p + 2X.r - 2 ^ CO) dz~\ (0 ^ .r ^ J g 1) (59) 
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where 

h 4 (x) = 

ti(x) > (0 £ x <£ £ 1) 

and X equals the largest positive value for which the quantity in the 
above brackets has one zero in x in the interval ^ x £ 1. Such a value 

must exist since / C(z)dz (x g x £ 1) is monotonic increasing. In 

addition, for this zero condition to hold independently of k, X -» /c(2p) 
and .€• — > 1 as /c — > ; thus, 



h.4(x) 



|~2 P - 2 f C(z)dz~\ (0 ^ x ^ x). 



(60) 



Test function fo(x) is seen to satisfy the same requirements as those 
discussed in connection with (52) through (55). Accordingly, interior 
bounds ±& are given by the relation 

£ < £ 3 = inf G(x,h) = inf (x + ft £j + 2X.r - 2 / C(z)<fe \ 

X<* I _'' : " J ° J J (61) 

- X-*Jb(2p)' = h 

fc-»0 

which is also consistent with the initial condition 

MO) = \^ 9(0) = |. 

Inasmuch as ft ^ £+ ^ £2 and & -» &(fc -» °)> the results of ( 58 )> ( 60 )> 
and (61) yield 



^ = |^A:(2 P ) ! (*->0) 



where 



f/0 (.r) ~ [2p - 2 J[ C(«) rf«J (ft 
P - f C(z) <fe. 



(62) 



0) 



Therefore, the choice of asymptotically equal test functions results in 
coalescent bounds. 
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In considering large values of A - , we take 

h(x) - x + X - x (O^igHKl) (63) 

where 

^ C(x) < 1 (OgK .£ ^ 1) 

C(x + X) -> 1 

X-.0 

x > a; + X > 

and X < according as (dC/dx)z I 0. Clearly, X is employed only in 
the event that C(x) has a discontinuity at x = x. Since the values of 
£ specified by 

£ < G(*A) = k ^ + 2tti + C(x) 
ax 

also satisfy the basic initial criterion 

kg(0) = £ < G(fiM0)) = -{x + X) + fcfc(O) ^ fcfo(0) 

interior bounds ±£4 are given by 

£ < £ 4 = inf <7(.r,/i 5 ) = inf [(fc - l)(x + X - *) + C{x)\. (G4) 

Z<£+X 

Employing the limits 

G(x,h b ) -» x (.r < 4 4- X) 

G[(x + X), fe(4 + X)] = C(x + X) 
yields 

Lim Lim £ 4 = Lim C'( i; H- X) = C(x) = 1 ^ £+ . 

X-.0 *-» X-0 

However, according to (23), £ + ^ 1; therefore, 

$+ = |_->l. (65) 

One can show in a similar manner that if the system captures for a 
given £ and k, a capture condition results for all smaller values of A\ 
Consequently, the relative capture as a function of A --2 assumes the 
general form depicted in Fig. 5 with the asymptotic behavior given by 
(02) and (Go). 
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Fig. 5 — General form of relative capture as a function of k 2 . 
VI. SAWTOOTH COMPARATOR — RC FILTER 

Directly related to the formulation and results of Section V is the RC 
filter-sawtooth comparator system in which 



C(x) = £ C (x -2n) (-co < x < oo) 
C (x) = 



n— * 

fx(-i <x < i; 

lo I.I-I ^ i 



(66) 



fm - T, fc = I — 1 , 



0"! 



= 1. 



The relatively simple form assumed by C(x) in the first domain en- 
ables one to derive exact expressions for both g (x) and £ . It is shown 
in Appendix C (cf. (99)) that 

y=J (k < 2) 

")\ (67) 



e. 



& - €f - I = 



k + exp 

1 (A- ^ 2) 



— ( tan ' T - 
to \ K 



where r 2 = (4 - A' 2 ) and < tan ' (•) < ir/2. For small A 1 , (67) can 
be expressed as 



ffl 



6. ~ * = 2 v i \*'i^° 



0"! 



(68) 



which checks with (62). In addition, from (20 et seq.) 

y P £ & . (69) 

An exact derivation of y p valid for the sawtooth and general filter case 
has been obtained by A. J. Goldstein. 6 
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As regards steady state phase error, relations (21) and (50) give 

X* = tT'ttdlwsi = |. (70) 

VII. SAWTOOTH COMPARATOR — LAG FILTER 

The capture results as developed in Section V do not apply com- 
pletely to the lag filter; however, the calculation of exterior and in- 
terior bounds in this case follows closely that of (40), (41) and (52). 
For B(/) = f, > 1, and C(x) specified by (66), (14) reduces to 

£ = C(x) + 0kg + g (l j g - - 2(0 - \)kg £ 8\x - (2n + 1)} 

dx V ( 71 ) 

- Qi(*a) 

where 
and 

fl-i + i. 

(To 

Using test function hi(x) of (38) yields the following noncapture values 
of£(cf. (40)): 

£ > Oi(xM(x ^ 0) 

where, as required, fcflf(O) = £ > fcj8fti(0) > khi(0). Hence, exterior 
bounds ±£ r are given by 

£,. = sup Gi(x,hi) = sup Gi(.r,Mi) 

£(4-/*V)« (M<V2) ( 72 ) 

1 (0* ^ \/2). 

Using separatrix solution 0/,(.'O of Appendix C as a test function yields 
the following capture values of £ (cf. (52) and (98) et seq.): 

£ < (-n(.v,y L ) = U (0 ^ x ^ 1) (73) 
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where 

7T-F =Tj (jSfc < 2) 

u - * + & - *' " Dl 1 -p (^ _ ta "^TJi - 'J/ 

1 (jSfc ^ 2) 

Tl 2 = (4 - /H- 2 ) 
and 

< tan _1 (-) < tv. 

Therefore, 

b^l^fe. < 74 > 

With regard to small fc, the relations for ^ and f L become 

fc,~/3fc (fe-»0) (75) 

fc~A; (fc-*0). (76) 

As reflected by parameter 0, the spread between £ v and £ L for small fc 
appears significant, but in most cases involving small k, a t » a and 
/3 = 1. Under more general conditions these bounds can be refined by 
constructing more complex test functions. 

According to the discussion in Section II (cf. (20) et seq.), £ f ; serves 
also as an upper bound on the pull-in range; viz., 

7p = 27+^2^. (77) 

VIII. SINE COMPARATOR — RC FILTER 

The sine-RC system is perhaps the most important practically, having 
received a great deal of attention in the literature. Several investigators 
have obtained capture results based on phase plane constructions of 
g (x) for large and medium values of fc. 1 J For small k, approximate re- 
sults have been derived by both analytic and experimental means, al- 
though these approximations are found to differ somewhat. 1 - 2 In this 
section we discuss the exact asymptotic form of the relative capture 
range for small fc. 

The normalized sine comparator is represented by 

C(x) = sin (irx) (-00 < x < oo ) (78) 
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whence 

p = f C(x) dx = - 

J IT 

/- = 1 

and 

Consequently, the first relation of (62) can be written as 

^f-^v': 1, ('••^°) (79) 

or 

^l/f,- 2 fe- )- «» 

The numerical limit in (80) has been estimated by McAleer 2 and Je- 
lonek 1 to be 1.00 and 1.27, respectively. As lower bounds, these approxi- 
mations differ appreciably. 

IX. SUMMARY 

The more significant results of the above development can be outlined 
as follows: 

(i) The variational formulation given by (27), a generalization of 
which allows one to treat similar problems in re-dimensional phase space. 

(ii) Bounds given by (42) and (61), yielding respectively noncapture 
and capture ranges for the general symmetric comparator and the sim- 
ple RC filter. 

(Hi) The first domain phenomenon associated with the symmetric-RC 
case of (ii) (cf. (56) et seq.). 

(iv) The exact asymptotic capture relations given by (62) for the sym- 
mctric-RC case. 

(v) The pull-in criteria given by (20 et seq.), (42), (61), and (62) for 
the symmetric-lag and RC cases. 

(vi) The specific capture relations given by Sections VI through VIII 
for the sawtooth-sine and lag-RC cases. 

(vii) The asymptotic relation given by (80) for small k, simple RC 
filters, and sine comparators. 
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APPENDIX A 

Critical Point Type 

To determine the behavior of trajectories near a critical point 

(x = x , g = 0), 

we linearize (11) and (12) about this point and investigate the natural 
modes of the corresponding linear solutions: As shown by Liapounoff, 5 
these modes identify the type of critical point in that the presence of at 
least one exponentially increasing component indicates a saddle point, 
and the presence of only decreasing components, a stable (nodal or focal) 
point. Equations (11) and (12) are linearized by expanding the right- 
hand sides in a double Taylor series in x and g about the point 

(x = .to, g = 0), 

and dropping all terms beyond the first order. Accordingly, we obtain 
from (11), (12), (21) and (22) 



(81) 



where 





dx _ 
dt 


ag 






dg _ 
dt 


bix 


+ big 


a = 


k 






hi = 


<T X (dC\ 

k Wxo 




b 2 = 


Tl[(l " 


-&) 


(a.-] 


ibine 


to yield 








d 2 x 
dt 2 


, dx 
bi dt 


— ab\X 



(82) 

the exponential components of the related solutions increasing or de- 
creasing according as the roots of the secular equation 

X 2 - b 2 \ - obi = (83) 
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are positive or negative. Moreover, if (dC/dx) Xo > 0, then 61 < 0, 
hi < 0, and 

X = !| 2 ± \ h j- + ahl < 0. (84) 



Also, if dC/dx Zu < 0, then 61 > and at least one value of X given by 
(84) is positive; i. e., 



x = "i + [£ + *J > 0. 



Hence, the critical point in question is stable provided (dC/dx) Xo > 
and a saddle provided (dC/dx) Ilt < 0. 

APPENDIX B 

Comparison Theorem 

Theorem 

If for a branch trajectory g(x(t)) and a test trajectory h(x), 

^ = F(x,g) (x eIJeJ) 

'£< {>)F(x,h) (xel) 

h(x„) ^ (^)g(x a ) 

h(x), g(x), t(x)eC (xeI,teJ), (85) 

where / denotes a closed interval of x(x„ ^ x ^ x b ) and J, a branch 
interval of t defined by the inverse function t = t(x) (x e I), then 

h(x) < (>)g(x) (x a < x ^ x„). (86) 

Proof 

Assume first that h(x a ) < g(x a ) and dh/dx < F(x,h) (x c I). With 
// and g continuous, either h(x) < g(x)(x e I) or the two trajectories 
intersect at some point .f in the interval, i.e., h(x) = g(x)(x a < x 
^ Xb)', however, an intersection implies that (dh/dx) x ^ (dg/dx) z , a 
condition which contradicts the combined relation 

©.<«*(»»-««*«)-(£).. (87) 
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Therefore, in this case, 



II 



and 



h(x) <g(x) (xel). 
h(x a ) = g(x a ) 



^<F(x,h) (xel) 

ax 

Consequently, h(x) < g(x) over some portion of /. Ruling out addi- 
tional intersections by (87), we obtain the general result 

k(x) < g(x) (x a < x ^ x b ). 

The alternative set of inequalities is shown to be valid in a similar 
manner. 



appendix c 

First Domain Trajectories in the Sawtooth Case 

For treating the sawtooth comparator-lag filter case, we seek here 
solution trajectories of (71) which span the initial state line and saddle 
points in the first domain, viz., %. x <: 1. As defined in (6(>), 

\x (0 ^ .i- < 1 ) 
C(x) = (88) 

[0 (x = 1). 

Hence, by (21 et seq.) the only saddle present in the first domain is that 
at the point (x = 1, g = 0). We intend, therefore, to determine a solu- 
tion g L (x) spanning the two points (0,£/fc) and (1,0). It is noted that 
in the simple RC case, g L (x) is identical to g (x) and £,. , to & (cf. (43) 
and preceding). 

Restricted to the interval (x ^ x ^ 1) and solution g L (x), (71) 
becomes 

iL = x + pkg L + g L ( ^-2({i-l)lcg l .5(x-l) (O^^l). (89) 
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However, since Ql(x) is required to vanish as x — > 1, gi,8(x — 1) = 0; 
consequently, 



h = x + 0ty L + g, ^ . 

a.r- 

Setting z = x — £/, and v = 0^/2 reduces (90) to the form 

2v dv 



(90) 



2— = - 

2 " i; 2 + 0kv + 1 



2v + |Sfc 



[91) 



FW 



... , , n_ y 2 + /3A'y + 1 

■ J '"4 



where 



n" = (4 - A 2 ). 
Integrating, we have 

In for*, 2 + fik(x - h)gt. + U - h)*\ 

~2g L + pk(x - h)' 



2/3/,: _, 

- 1 — tan 



Ti(.r - fc) 



+ C 3 (92) 



where 0» = const. At the saddle point (g L = 0, x — 1), (92) becomes 
|n(l -^-^tan-'^ + C.. (93) 

Therefore, <7 t is given implicitly by the relation 

r/, 2 + pk(x - u)9l + (.1- - h) s 



= cxp 



(1 - fc) 1 

2/3/,- 



7"! 



, -, /2f/, + jSfc(z - £,)\ _, //3/,- 

tan I — I - tan 

Tl(.V ~ &,) 



*"l /_ 



(94) 



An explicit expression for £/. is derived by inserting the initial condition 
{//.(0) = fc/fc in (94); i. e., 



h:\ \ -):-{&- 1) 1 
/.--'L (1-fc) 1 J 



-{?[^(^) --"'©]}■ 



(95) 
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For £ L ^ 1 (cf. 23) and 0k < 2 (cf. 91 ), the argument of the first arc-tan 
function in (94) takes on an infinite value at only one point in the 
interval (0 i x ^ l);i. e., 



2g L + 0k(x - g t ) 
ri(.r - h) 



fik 



> 



— oo 

/3fc 2 - 



(s-1) 

(x -> h + ) 

(x -> tn 

tj^)jo (—0). 

Thus, the arc-tan difference in (95) must be such that 

-<-i^) —-©--- '[*=%*]<» <»> 

Equation (95) can then be written as 

£J1 -^-^=-{f[--(^v t )]} (9 



A'(l - &) 



7) 



or 



«L = 



* + [1 - tfifi - Dl 5 exp g' [tan" 1 (jg^jj) ~ *]} (98) 

(jBfc < 2) 

where < tan _1 (0 < *• Thatfc, = 1 forjSfc ^ 2 is shown by considering 
two values of file, say (Sfei -» 2 and 0fa > 2, and two trajectories with re- 
spective initial conditions fi/fci = 1/A'i and | 2 /A- 2 = 1/Ar 2 . Since from 
(98), gt — » 1 = 6 = & as jSAi — > 2, the former trajectory is a capture 
type. However, 1/fci > l//c 2 ; therefore, the latter initial point and 
trajectory lie below those of the former. This implies that relations 
£ 2 = 1 and /3/c 2 > 2 correspond to capture conditions. On the other hand, 
if & > 1, noncapture trajectories result. Hence, £ l = fe = 1 for /3fc > 2. 
Finally, capture criteria in the simple RC case are given by 



£t "TF*r* £o — 



A- 



fc + exp 

1 



j=— v*i (A^ < 2) 

(fe ^ 2) 
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where 



and 



APPENDIX D 



to 2 = (4 - tf) 



< tan l ( • ) < tt/2. 



List of Symbols 

(Pi — input phase 
ipa — output phase 
e — <pi — <p — phase error 



x = normalized phase error 



TV 



Xtt — normalized steady-state phase error in the syn- 
chronous state 

/(e) — comparator output 

/ m — sup | /(e) | 

C(x) = ' — normalized comparator function 

J in 

P = f C(z) dz 

•Ml 

u c — free-running frequency of the local oscillator 
Aco(0 — input frequency deviation 

Acoit) 
Sl(t) — ~r~ — normalized input frequency deviation 

&J in 

C\ = const. — magnitude of input frequency step 

Ci ... 

£ = . — relative input frequency step 

£+ , £- — relative positive and negative capture ranges, re- 
spectively 
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£ r = £+ — £_ — relative capture range 

go — relative input step corresponding to separatrix 

go(-v) 

C< = const. — magnitude of asymptotic input frequency step 

y = — r — relative asymptotic input frequency step 
afm 

y+ , y _ — relative positive and negative pull-in ranges, re- 
spectively 

7( , = 7 + — 7_ — relative pull-in range 

H(t) — filter impulse response 

S(t) — filter output 

= h_ dx 

9 aidt 
g (x) — primary first domain separatrix, as defined by the 
description prior to (43) 

5(t) — Dirac delta function 

a — gain constant; viz., — = co, : + aS(t) 



k = 



IT (J i 
5 



<p! , c-i — filter parameters (cf. Fig. 1(b)) 
— normalized system parameters 



= 1+^= 

* — convolution 

s — Laplace transform variable; i.e., £•/(<) = F(s) 
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